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^ ■ ABSTRACT 

a: 



This series of papers investigates the dynamic interior of a quiescent prominence 
revealed by recent Hinode and SDO/AIA high-resolution observations. This first 
paper is a study of the static equilibrium of the Kippenhahn-Schliiter diffuse 
plasma slab, suspended vertically in a bowed magnetic field, under the frozen- 
C/} ■ in condition and subject to a theoretical thermal balance among an optically- 

thin radiation, heating, and field-aligned thermal conduction. The everywhere- 
analytical solutions to this nonlinear problem are an extremely restricted subset 



of the physically admissible states of the system. For most values of the total mass 
frozen into a given bowed field, force-balance and steady energy-transport cannot 
both be met without a finite fraction of the total mass having collapsed into a 
cold sheet of zero thickness, within which the frozen-in condition must break 
down. An exact, resistive hydromagnetic extension of the Kippenhahn-Schliiter 



> 

slab is also presented, resolving the mass-sheet singularity into a finite-thickness 
layer of steadily-falling dense fluid. Our hydromagnetic result suggests that the 
narrow, vertical prominence H a threads may be falling across magnetic fields, 
with optically-thick cores much denser and ionized to much lower degrees than 
conventionally considered. This implication is discussed in relation to (i) the 
recent SDO/AIA observations of quiescent prominences that are massive and yet 
draining mass everywhere in their interiors, (ii) the canonical range of 5 — 60 G 
determined from spectral-polarimetric observations of prominence magnetic fields 
over the years and (iii) the need for a more realistic multi-fluid treatment. 

Subject headings: MHD -- Sun: magnetic fields -- Sun: corona — Sun: promi- 
nence 
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Introduction 



Quiescent prominences are long-lived, part ially-ionized, plasma condensations in the 
tenuous, million-degree hot, ionized solar corona (jMartin. Bilimoria fc P. W. Tracadad Il994l ; 



Tandberg-Hanssenl Il995l ; iGaizauskas! Il998l ; lLabrosse et al.l I201CH : iMackay et al.l 120101 ) . These 



partially- ionized condensations, at least, their optically-thin observable parts, are typically 
two orders of magnitude denser and cooler than the surrounding corona. As a macroscopic 
structure prior to an eruption, a prominence can be remarkably stable, persisting for days 
to weeks as a long, vertical slab, approximately 5 x 10 3 km wide, 2 x 10 4 km tall, 10 5 km 
long, seemingly suspended by magnetic fields to curve the full prominence length right above 
and along a magnetic neutral line on the photosphere. By prominence we shall henceforth 
be referring to this common type, the polar crown prominence being the canonical example; 
see the above review articles for a comprehensive view of prominences of all varieties. 

Figure 1 is a Hinode/SOT image taken in the Ca II 396.8 nm "H-line" spectral region 
showing the cool part of the prominence at a high spatial resolution, revealing its characteris- 
tic narrow, vertical, dense threads interspersed among narrow, dark lanes. This static image 
belies the remarkably dynamic nature of these interior structures. Prominence movies from 
Hinode and SDO/AIA instruments have established that over internal scales of the order 
of 300 km and time scales of minutes, these structures are in a constant state of motion at 



less than free-fall speeds in a characteristic range of 5 — 30 km s 1 ( JBerger et al.ll2008l 12010 



2011uOkamoto et al.ll2007t lOkamoto. Tsuneta fc Bergerll2010l ). The ever-present cool, dense 
threads are falling whereas the dark lanes amon g them are narrow strea ms of hot, tenuous 
plasma traveling upward at comparable speeds (JBerger et al.ll2008U2010l ). Vertical motions 
in both directions are actually present in both types of structures, but the high-resolution 
movies show a discernible average descent of the dense threads and a general bubbly rise of 
hot, tenuous plasmas. The latter may be the result of rapid heating to coronal temperatures 
associated with emerging magnetic fluxes that push into the overlying dense prominence 
material. This condit ion may be understood in terms of a n onlinear Rayleigh- Taylor hydro- 
magnetic instability ( jBerger et al.ll2010t iHillier et al.l 120111 ) . Such emergent-flux structures 
can be of considerable sizes that first form quasi-statically beneath a quiescent prominence 
before pushing their ways clear through the prominence without a macroscopic eruption 
rtde Toma et~aDl2008l lierger et alJboilh . 



A recent quantitative analysis of the motions of descending threads observed with 
SDO/AIA indicates that over a 20 hour period, an order of magnitude more mass may 
be drained i n this manner t han found at any one time in a moderate-size, slowly evolving 
prominence ( JLiu et al.ll2012l ). Using the lower end of the range of observed prominence den- 
sities to compute, the drainage can remove some 10 15 g from the corona in a day, a mass 
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comparable to that of a coronal mass ejection. Of particular interest is that the drainage 
rate estimated from observation is changing in proportion to the slowly-changing, estimated 
total prominence mass during quasi-steady evolution. This observational result raises two 
questions. How does the observed significant mass drainage via cool threads relate to the 
replenishing mass-injection represented by the upward streams of hot, tenuous plasmas? 
Berger et al. (2011) propose that emerging magnetic flux drives a constant cycling of mass 
from the photosphere, in the form of tenuous hot plasma, into the coronal cavity overlying 
a quiescent prominence. This up-flow returns as the draining vertical threads, essentially 
establishing a novel form of 'magneto-thermal' convection in the solar outer atmosphere. 

The other question relates to the embedded magnetic fields in these processes. Spec- 
tral polarimetric observations have shown that quiescent prominence plasma s embed typ- 
ically horizontal, local fields of the order of 5 — 60 G (ICasini et al.l 120031 : iLeroyl 11989 



Lopez Ariste fc Casinill2003l ; lTrujillo Bueno et al.ll2002l ). These fields are sufficiently intense 
to dictate the structures of their prominences under a high degree of the frozen-in condition 
due to high electrical conductivity. In particular, a large amount of mass can be supported 
by fields of the observed intensities. Yet more mass can drained through the prominence 
in less than a day than the amount trapped quasi-steadily in the prominence at any one 
time. How is the observed drainage reconcilable with the observed prominence field being 
generally horizontal if we accept the frozen-in condition? 

An attractive hypothesis is that it is in the nature of prominence dynamics to develop 
such extreme physical conditions that the usual frozen-in condition breaks down sponta- 
neously. This hypothesis is what we set out to investigate as a theoretical hydromagnetic 
possibility. The first paper of the series here treats the case of thermal collapse during 
plasma condensation producing cold, dense, poorly ionized plasmas within which the frozen- 
in condition cannot be sustained. The second paper takes up a separate piece of physics, 
demonstrating that despite the high magnetic Reynolds-numbers of the prominence inte- 
rior, magnetic reconnec tion takes place readily via spontaneous formation and dissip ation of 
electric current- sheets (jParkerl Il994j ; IPetrie &: Lowll2005l ; IJanse. Low fc Parker! I201CH ) . The 



two papers provide a hydromagnetic basis to understand the observed restless interior of a 
prominence. 

Specifically the first paper here treats the following theoretical problem. We construct a 
complete family of exact solutions describing a vertical Kippenhahn-Schluter (KS) plasma- 
slab supported against gravity in a bow ed magnetic field in unbounded space, varying only in 



the d irection perpendicular to the slab (IKippenhahn fc Schuterlll957l ; lLowlll975l ; iHillier et al. 



201ll ). We assume static force equilibrium under (i) the frozen-in condition and (ii) a steady 
balance among optically-thin radiation, heating, and thermal conduction directed along the 
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field. These three energy processes take the theoretical forms studied by Low & Wu (1981), 
representing an energy transport of a minimum physical complexity relevant to prominences. 
Despite its one-dimensionality, the KS slab captures important properties of this complex 
system. The novel result of the study is that the set of everywhere-analytical solutions 
found in the earlier study is only a subset of measure zero of all the physically admissible 
equilibrium states of the KS slab. Generally, force-balance and steady energy-transport 
cannot both be met under the frozen-in condition without a fraction of the total mass of 
the slab having collapsed into a mathematical singularity representing an exceedingly cold, 
dense sheet at the center of the otherwise diffuse slab. The weight of that mass sheet is 
supported by a discrete Lorentz force due to a discrete electric current in the sheet. This 
discrete current can be maintained only if the frozen-in condition is rigorously obeyed, that 
is, the electrical conductivity is infinite. The electrical conductivity in the solar atmosphere 
is large but finite, of course. In that case, an evolution to equilibrium must at first seek 
a state with an inevitable singularity when the frozen-in condition does apply to a very 
high degree. As the singularity develops, the small nonzero electrical resistivity unavoidably 
becomes significant at a point in time, resulting in current dissipation. A breakdown of the 
frozen-in condition thus occurs spontaneously. 

The general inevitability of mass-sheet singularity in the KS-slab is demonstrated in 
Section 2. The singular equilibrium solutions containing a mass sheet with a discrete current 
are analyzed in Section 3. The breakdown of the frozen-in condition is illustrated in Section 
4, treating an exact solution for a steady, vertical resistive hydromagnetic flow and pointing 
out the need for a more realistic multi-fluid description. Our results are summarized in 
Section 5 to relate in general terms to the physics of prominences . 



2. Coupling between static force-balance and steady energy-transport 

We adopt a one-fluid hydromagnetic description of the corona and prominence: 

= i(VxB)xB-Vj)-^, (i; 

An 



dv 


■V)v 




dB 




~dt 



: Vx(vxB-j|VxB), (2) 

coupling the velocity v to the magnetic field B; p, p, and g being the density, pressure and 
uniform gravitational acceleration in the negative Cartesian- z direction. The coefficient of 
resistivity rj = c 2 /Alter is generally very small over all observable length scales in the sense 
of a large magnetic Reynolds-number Rm = VqIq/t), c and a being the speed of light and 
electrical conductivity, respectively, in cgs units, and, vq and Iq denoting characteristic speed 
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and length, respectively. The time To = Iq/tj for a magnetic structure of characteristic 
length Iq to diffuse resistively is another indicator of the importance of resistive effects. Set 
l m 7 x 10 7 cm as the limit of spatial resolutions achievable with a ground-based telescope. 
Then the Spitzer (1962) resistivity r\ = 3 x 10 12 T~ 3 / 2 cm 2 s^ 1 , where T is the temperature, 
gives To well in excess of tens of years or longer for temperatures in the range of 10 4 — 10 6 °K. 
Therefore, the frozen-in condition can only break down if a hydromagnetic fluid develops 
exceedingly small scale structures. In the case of the corona, these scales may even fall 
below the mean free path of the fluid, so that the app licable resistivity involves coupling 
between the large-scale fluid and particle- kinetic effects ( ILow et al.ll2012l ). 



We are interested in a static equilibrium with a negligible r], so that currents may persist 
permanently. Then, setting v = 0, the induction equation (J2J) is irrelevant except to demand 
that B is solenoidal. The mass- conservation equation is also irrelevant, and static equilibrium 
is described by the complete set of equations: 



k- 



4tt 
(B ■ VT) B 



(V x B) x B - Vp - pgz = 0, 
V ■ B = 0, 



IBI 



-r(p,T) + h(p,T) = 0, 

V = — pT, 

m 



(3) 
(4) 

(5) 
(6) 



where we introduce the ideal gas law, with ks and m being the Boltzmann constant and 
the mean mass of the fluid particles, respectively. Force-balance equation flH]) is coupled to 
steady energy equation (J3J. The latter balances radiative loss r, heat gain h and a thermal 
conductive flux aligned along B with a thermal conductivity k. The properties of a specific 
plasma fluid is represented by r, h and k as known functions of the thermodynamic variables, 
p, p and temperature T. 

This radiative magnetostatic problem is formidable, especially for a three-dimensional 
system. We therefore treat a localized one-dimensional system. To keep physical issues in 
mind we first briefly describe the model in relation to the real prominence and its surrounding. 
Then we formulate the mathematical problem in Section 2.1 and show in Section 2.2 that 
the physically complete set of solutions must include both everywhere-analytic solutions as 
well as the solutions with admissible singularities presented in Section 3. 

The prominence in Figure 1 belongs to a larger-scale, three-part coronal helmet-streamer 
of the kind in the white-light image of a total eclipse in Figure 2. Coronal o bservation of the 
magnetic fields external to prom i nences in the corona is still in its infan cy ( jCasini fc Judge 
19991 ; lJudge. Low fc Casinill2006l ; iTomczyk et al.ll2008t iDove et al.ll201ll ). but it is physically 
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reasonable to interpret a helmet-streamer to be an arcade of closed magnetic field anchored 
with magnetic footpoints frozen into the base of the corona (JLowlll996l . l200ll ; IZhang &: Low 
20051 ). This field traps the high density of the helmet-streamer in approximate static equi- 
librium under the frozen-in condition. External to the helmet-streamer is the solar wind 
accelerated from quasi-static conditions at the coronal base. This wind combs its frozen-in 
fields into an open configuration, each field line extending from one footpoint anchored in 
the coronal base out into the heliosphere. The bright helmet often has a tenuous, dark cavity 
at its base within which a quiescent prominence can be found located at the bottom of that 
cavity. When the line of sight is aligned favorably with the arcade of the helmet-streamer, 
an exemplary case in Figure 2, this three-part structure is seen clearly. 

Observations support the interpretatio n that the cavity is a magnetic flux rope with field 
lines that wind around with in that cavity. ( jFan Il200ll ; Manchester et al.ll2004J : IZhang fc Low 
20051 ; ICottaar fc Fan! 120091 ) . The prominence may then be explained as the cool condensed 
plasma trapp ed at the bottom of l ocal U-shaped flux tubes under a high degree of the frozen 



in condition (IFan fc Gibson 



2006 



200 



7; 



Fu 



ler et al. 2008; Gibson fc Fan 2006; Gibson 2010 
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Schmit et al.ll2009t iSchmit fc Gibson! 120111 ). The frozen-in condition permits no motion of 



the condensation across magnetic flux surfaces. Our study here sets out to first understand 
this frozen-in state well, only to discover from the physics of such a state that the frozen- 
in condition must generally break down in the prominence interior to produce cross-field 
transport of mass at less than free-fall speeds. 



2.1. A thermally-balanced Kippenhahn-Schluter slab 



Consider a ID system that depends on only the Cartesian coordinate y. Let the field 
have the solenoidal form: 

B = B [Go,l,H(y)], (7) 

Bq and Go being constants. We are interested in B z = BoH(y) being a monotonically- 
increasing, odd function with H = at y — 0. This field is deformed from a uniform 
"background" field 

Bfc = B [G? ,l,0], (8) 

by weighing upon it with frozen-in mass, producing a bow-shape geometry to represent a local 
U-shaped. In this idealized unbounded domain, the surrounding prominence environment 
tends to a uniform state located at y — ¥ ±oo. 



Introduce the dimensionless variables, 



P = p/p , D = p/ Po , 9 = T/T , Y = ^- 



(9) 
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where Lq = -^-^ is the hydrostatic density scale height at temperature T related to the 



rung 



other normalization constants by po = pokTo/rrio. The ideal gas law is then 

P = DO. (10) 

Introduce the two plasma /9-parameters: 



Bl 



P = ^ (") 



87Tgo 

Pnet ~ 5 2 (1 + H{Y) 2 + Gg) P ' l j 

respectively measuring the importance of the fluid pressure relative to the magnetic pressures 

1 /9 

of the uniform B y = B and the total field |B| = B (1 + H 2 + Gfy . With no loss of 
generality, take p = po to be the pressure at the center of the sheet, ie, P — 1, H — at 
F = 0. Then, the force-balance equations ([3]) reduces to the y and z components: 

(3P + H 2 = (3, (13) 

§->■ 

Equation (TTB^) shows that the total pressure p + B 2 / 8tt is uni form in space which is a sufficient 



condition for linear hydromagnetic stability (jZweibelllf982l ). The high pressure of the sheet 
is confined by the vertical component B z = B H significant outside the sheet. Equation (J3J) 
can be written in the form 

-L(B-V)B-v(p + !^ -pgz = 0, (15) 

in terms of the magnetic tension and the total-pressure force. Since the total pressure of the 
KS-slab is uniform in space, the tension force is everywhere vertical and it alone supports 
the weight of the plasma. In other words, the left hand side of Equation (FH|) is essentially 
the upward magnetic tension force. 

We pause to relate the above idealized magnetic field to spectral-polarimetric obser- 
vations of prominences. If we identify the concentration of plasma around Y = with 
the prominence plasma being observed, where B z = B H w 0, polarimetric observation 
will detect the horizontal field B& fl given by Equation (jSJ). The observed field has intensity 
Bq (1 + Gq) which can fall within the observed range of 5 — 60 G even with Bq small as 
long as Go is sufficiently large. Therefore, the plasma-beta /3 given by Equation (fTTjl may be 
larger than unity whereas the plasma-beta (3 net defined by Equation (Tl2l) in terms of the net 
field is typically less than unity. This feature simply reflects the fact that a slab of plasma 



hemmed up by a background field aligned close to the central plane of the slab will snag 
only a weak component B y of the field across it. Therefore, the field would sag greatly in 
the y — z plane with a large value of j3. 

We take advantage of this simplicity of the model to deal with the complexity of thermal 
balance. We adopt the temperature-dependent Spitzer (1962) thermal conductivity 



Ks 



rp5/2 



erg cm 



■i K - 



(16) 



for a fully ionized gas where kj = 2 x 10 6 . For our theoretical purpose, we adopt this 
conduction model in a formal sense with ks taken as an ad hoc constant, so as to be 
conceptually specific about the model of thermal conduction adopted. To complete the 
specification of the fluid's material properties, we take 



r 

h 



a p 2 T, 
7oP, 



(17) 
(18) 



where ao and 70 are constants characteristic of the fluid. Here we model radiative loss 
r to be optically thin, collisional in nature as represented by its quadratic dependence on 
density. Detailed radiative calculations have shown that the collisional r has this form with a 
approximately constant in contiguous ranges of temperatures from 10 4 if to 10 6 K, multiplie d 
by a monomial in T of a constant power in each of those ranges ( JTandberg-Hanssenl Il995l ) . 
In Equation (TT7|) we take a® to be constant and the temperature factor to be linear in T 
to simplify an otherwise very complex mathematical problem. The origin of heating in the 
prominence is poorly understood, representing another uncertainty in such th ermal models . 
For our theoretical study we use an ad hoc heating rate proportional to density rtUchidalll963l: 
Tandberg-Hanssenlll995f ) . These choices lead to an analytically integrable model (ILow &: Wu 
198ll ). We wish to investigate this model afresh to understand the physical meaning of these 
analytical solutions. 

With the material properties of the plasma fluid so specified, energy equation (J5J) takes 
the dimensionless form: 



d 

dY 



K e 5 / 2 dd 

_l + G 2 + H 2 dY 



D 2 6 - iD, 



where 



7 



Kq 



70&B 87T7ofcs 1 



2 ft' 



a m p a m B^ f3 



KsTc 



5/2 



2 r2' 



«OPo L 



(19) 

(20) 
(21) 
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L being the scale height at temperature T already introduced. Equation ([19]) as a second 
order ODE requires two boundary conditions to determine 9. We choose them to be 

d6 7/2 
Y -0, ^ = 0, (22) 

d9 7 ' 2 , s 

Condition (|22|) comes from the assumed symmetry of 9(Y) about Y = 0. Since the heat- 
ing source is explicitly accounted for by h, condition (T23"j) sets the thermal flux to vanish 
at infinity. In the unbounded model considered here, the /i-heating and r-radiation vanish 
at infinity. An external source of thermal flux coming in from infinity would require an 
unphysical temperature that increases to unbounded values at infinity. An external sink of 
thermal flux at infinity is incompatible with the unbounded domain, since a monotonically 
decreasing temperature must become zero at some finite Y on either sides of Y — 0. The 
/i-heating is the only source of energy in this magnetostatic atmosphere with thermal con- 
duction transporting excess heat to where it can be radiated away. The equations (fTOl) . (TT3|) . 
(TT4|) . and (TT9l) are a complete set of scalar equations for the dependent variable P, D, 9, and 
H . In principle, Equations (II 01) . ( 1T3|) . and (1141) determine P, D and H in terms of 9 so that 
Equation (I19p determines 9 as a function of Y subject to the two homogeneous boundary 
conditions ( E2|) and (123"j) . From the solution 9{Y) obtained, the other dependent variables 
can then be expressed as functions of Y. 

The model is ID only in the mathematical sense. The field is sheared and bowed in 
3D space. This 3D geometry strongly influences energy-balance, represented by the factor 
(1 + Go+i^ 2 ) -1 in Equation (1T9|) . The thermal flux is completely aligned along the field. This 
factor depends on the magnitudes of the two field components, B x = B Go and B z = B H. 
For large values of Go and H , respectively describing high magnetic shear and great sag in 
the field, the thermal flux enters the plasma slab at a large oblique angle. 



2.2. The physically complete set of solutions 

The governing Equations (!T0|) . (!T3|) and ( TH|) and ( TT9l) can be reduced to a quadrature. 
This was shown by Low & Wu (1981) using P as an independent variable and D, 9, H and 
Y the dependent variables. In that formulation, D, 9, H can be solved analytically explicitly 
in terms of P, following which solving Y as a function of P gives the distributions of the 
equilibrium field and plasma in space. Here we take a different approach, choosing H instead 
of P as the independent variable. 

By Equation ( 1T4"|) describing the vertical balance of force in a slab symmetrical about 
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Y — 0, H(Y) must be monotonically increasing in the domain — oo < Y < oo since density 
p = p D is positive. Define H (±00) = ±H , and Equation f|T3|) gives if — a/5 with 
P(±oo) = 0. The total mass M in an infinitely long column, of unit cross-sectional area, 
parallel to the Y-axis is given by 

/oo 
DdY 
-00 



^ 



2, • * (24 > 

Hereafter we called M simply the total mass. During an evolution of this ID system, M is a 
constant in time. When this system is in equilibrium in a background field with a given B , 
M determines /3 = 87rpo/-Bo ^y Equation (124p . In other words, M determines the central 
fluid pressure p(0) = po and the far vertical field component B z (±oo) = ±BqHo, recalling 
Hq = Vfi- The greater the total mass the greater is the sag of the field with a large far 
vertical component and a high pressure po at the bottom of the V-shaped field. 

We now pose the following mathematical problem. Equation ( fl3l) gives P{H) explicitly. 
If 9(H) is also known, D(H) is given by the ideal gas law (TTOl) . Then solving Equation (TT4l) 
for H{Y) not only determines B(F) but also gives P(Y), 0(Y), and D(Y). Therefore, the 
problem posed starts with determining 9(H) in the domain \H\ < Hq = y/]3 for a specified 
p. In other words, we need to transform energy equation ( 1191) as a differential equation for 
9 from domain \Y\ < 00 to an equation in the domain \H\ < Hq = -\//3. Let us carry out 
such a transformation and solve for 9(H) in the three-step development below, leading to 
a family of physically complete solutions, the meaning of "physically complete" becoming 
clear at the end of the construction. 



2.2.1. The isothermal KS slab and related simple models 

It is instructive to consider the isothermal solutions, treated by Kippenhahn and Schluter 
(1957). Such a state arises if the fluid does not radiate, cto = 0, and is not heated, 70 = 0, 
but is thermally conducting ks 7^ 0. Then with 9 = 9%, a constant, and P = 9±D, Equation 
(TJJJ can be integrated to give the isothermal solution: 

H = v^tanh^F, (25) 

D = ^ech 2 7^F. (26) 
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The field lines in 3D space has uniform B x = BqGq and each projects into the bowed curve 
in the Y — Z plane: 

Z - Z = 20! log cosh %£Y, (27) 

identified by its crossing the axis Y = at Z = Z in normalized unit L . This isothermal 
solution is displayed in Figure 3. 

Along the projected field line, the invariant frozen-in mass M is distributed hydrostati- 
cally, decreasing with Z exponentially with a scale height LqOi, recall that d\ is temperature 
in units of the free normalization constant T . This exponential decrease with height com- 
bined with the self-consistent shape of the field line renders the pressure P a strict function 
of Y, seen in the following two equivalent expressions, 

seen 2 ^ 

= exp(^), (28) 

the latter giving P varying with height Z along a fixed field line identified by the constant 
value Z . As a function of Y, P has a width proportional to the scale height L 6i and 
inversely proportional to the total mass M, with yfft related to M by Equation (liMl) . With 
the background field B& 9 fixed, the greater M is, the greater would be the sag in the field, 
and, hence, the narrower would be the peaks of P and its isothermal density D. 

On a scale large compared to the pressure width, we only see the two uniform-field 
parts of the bowed field lines on two sides of Y = 0. The mass concentration centered at 
Y = appears as a sheet of zero thickness on such a scale. On the other hand, letting Q\ 
go to zero decreases the scale height L$\ and the width of the mass concentration both to 
zero in absolute measure of length. This happens, for example, in a fluid that is not heated 
but constantly radiating away energy, with a very high thermal conductivity to keep the 
temporally decreasing 9\ uniform in space. If we fix /3 in this process, the total mass M is 
unchanged. In the limit d\ — )■ 0, D becomes proportional to S(Y), the Dirac delta function 
accounting for the unchanging mass M. At the same time, H given by Equation (1231) tends 
to a step function jumping from — y/]3 to y/]3 upon crossing Y = from the left. The pressure 
changes such that P = D9\ — )■ 1 at Y = for all time but vanishingly small for all points 
7^0, self-consistently describing a cold thin sheet of mass in pressure equilibrium with 
the external magnetic field. In these considerations physics dictates that fluid and magnetic 
pressures must always be bounded whereas density may grow unboundedly so long as its 
integral gives a bounded mass. This narrowing of the density peak with a fixed total mass 
combined with a vanishingly small temperature illustrates a form of thermal collapse relevant 
to Section 2.2.3. 
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Two related properties are worthy of mention when the temperature varies in space. 
Equilibrium requires B z = B H(Y) to be monotonically increasing in — oo < Y < oo. 
This ensures the pressure is positive and monotonically decreasing with Y. It also ensures 
the density is positive but does not require the density to be a monotonically decreasing 
function of Y. Depending on the variation of the temperature with Y, the density may 
have two maxima located off the center Y = 0, rather like the density inversion layer in 



a vertical atmosphere ( iLerche fc Lowlll980l ). The other property is that in the unbounded 
space, energy balance may be such that the pressure P declines with Y to become zero at 
some finite distance Y = ±Yq external to which P = D = 0. The mass sheet has condensed 
and contracted t o a finite width rath er like a vertical atmosphere with a finite top above 
which is vacuum ( ILerche fc Lowlll980[ ). 



2.2.2. The everywhere- analytical solutions of Low & Wu 

Let us transform energy equation (TT91) in F-space into its equivalent, 

(3-H 2 



10 °' dH 



de b ' 2 

l + Gl + H 2 dH _ 



0(l-7)-ff' 



(29) 



in if-space, using Equations (TTU1) . ( TT3l) and (fl4l) to remove explicit dependance on Y. The 
term within the square brackets on the left side of Equation (I29p is essentially the thermal 
flux along the field. Introduce the dimensionless thermal flux T defined by 



T = -(3 2 k 



£5/2 



d0 



l + G 2 + H 2 dY' 



(30) 



Removing the explicit Y-dependence from its right hand side, we rewrite this definition as 

d9 5 / 2 



J '~ 10 P K °l + G 2 + H 2 dH ' 



so that Equation ( 129|) takes the form 



with the integral 



g-fll-7)-*, 
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(31) 

(32) 
(33) 



Hi being an integration constant. A second integration using Equation (J3TJ) then gives 
3 
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(34) 



The homogenous boundary conditions (12211 and (1231) require J 7 = at H = and 
H = Ho = y/]3 in the if -domain, leading to H\ = and 7 = 2/3 and the solution 



where T = 1 + G^ + /3 and #2 is an integration constant. 
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the everywhere analytical solution of Low & Wu (1981). The free integration constant 62 
is the temperature at H = ±y/]3 corresponding to Y — > ±00. The temperature decreases 
inward monotonically to the temperature 9\ at Y = H = given by 

5 /.. ..-,> I 

2' 



6-> 



5/2 



eT - 



1 + G 2 + ^/3 



(36) 



3f3 2 K 

The integration constant ^2 is arbitrary but is physically required to be positive and bounded 
above such that Q\ > 0. Let us replace 82 with 9\ as the arbitrary integration constant, with 
no loss of generality. 

Our boundary value problem with its homogeneous boundary conditions is a nonlinear 
eigenvalue problem that admits only one value, 2/3, for the eigenvalue 7. We use the nota- 
tion 7 C = 2/3 whenever we mean 2/3 as a special number in this sense. The eigenfunction 

05/2 correS p 0nc [i n g to 7 = j c is highly degenerate, that is, there is an infinity of eigenfunc- 

5/2 
tions, corresponding to the second free integration constant 9^ not fixed by the boundary 

conditions. For each member of this infinite set of solutions in the ii-domain, defining the 

temperature 9, the governing Equations ( ITOfl . ( 1T3l) and ( 1L41) can be used to construct the 

corresponding spatial distribution of field and fluid. 

The eigenvalue 7 = 7 C is set by the global balance of energy. Homogenous boundary 
conditions ( 122]) and (123]) describe a closed thermally-conducting system. As we have pointed 
out, thermal conduction is not a source of energy in this unbounded domain, serving only to 
conduct heat to balance the total heat input 1-L against the total radiational loss 1Z. These 
two quantities can be explicitly evaluated in the if-domain: 



n 
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vr (37) 



with the ratio 
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/o 
-00 

- 3^' (38) 

I = ^ < 39 > 

The eigenvalue 7 = 7 C = 2/3 follows from "H = 1Z for global energy balance, described by 
a whole family of everywhere-analytical solutions distinguished by the temperature Q\ at 

Y = 0. 

Figures 4-6 show representative solutions for different values of /3 and Go, each with 
a temperature profile rising from some central temperature 6 = Q\ at Y = to a hundred 
times higher temperature = 82 at Y — > 00. This is rendered by picking a suitable value for 
the constant k,q characterizing the importance of the thermal conductivity to the radiative 
loss. The constants (/3,Gq) fix the amount of frozen-in mass M and the obliqueness of the 
field-aligned, thermal-conductive flux supplying energy to the central part of the slab. 

Figure 4 shows a solution with (3=1 and Go = for a field lying entirely in the 

Y — Z plane, a useful reference state for comparison with the equilibria in Figures 5 and 6. 
Stratification by gravity implies that radiative loss and heating, respectively, dominate in 
the dense region around Y = and in the far rarefied regions. The steady state is naturally 
one of high temperatures in these far regions to drive a pair of thermal fluxes inward to be 
radiated away. The narrowly peaked density at Y = is produced by this global balance of 
energy. The peaked density implies an enhanced radiative loss and a corresponding depletion 
of the incoming thermal- conductive fluxes. The temperature gradient is extremely steep near 

Y — 0, due in part to the power-law dependence of the thermal conductivity on temperature. 

Figure 5 may be interpreted as a solution with the same B with Go = as in Figure 
4, but with a greater amount of frozen-in mass M so that (3 = 5.0. The more pronounced 
sag of the more heavily loaded field is obvious from the projected field line described by 
Z(Y). The greater sag also implies a steeper rise in temperature as we leave Y = owing 
to the more oblique incidence of the thermal-conductive flux along the field. Alternatively, 
the strong B z = B H(Y) in this case may be said to be partial insulating the central part 
of the slab from the extremely hot exterior. 

The temperature gradient in the vicinity of Y = can be further enhanced by increasing 
the oblique angle of incidence of the thermal flux with Gq 7^ 0, illustrated by the solution in 
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Figure 6. We recall that /3 is the central plasma beta defined relative to the magnetic pressure 
due to the uniform component B y = B . If we keep the background field intensity |B& fl | = 
B (l + Gq) 1//2 fixed, we may increase G as B decreases and (3 increases correspondingly 
One way to think of this is to take the slab in Figure 4 of a fixed mass and juxtapose 
it to align the plane of the slab more with the background field of a fixed intensity. The 
solution in Figure 6 shows the case with the slab oriented such that B x /B y = G = 5, giving 
f3 = 130. The same slab mass as in Figure 4 is now supported in the deep sag in the field 
line projected on the Y — Z plane, as described by Z(Y). The increased insulation effect of 
increased tan~ 1 H and tan -1 Go gives rise to the extremely steep temperature gradients in 
the vicinity of Y = required for energy balance. 



2.2.3. The physically complete set of magneto static solutions 

The set of everywhere-analytical solutions are an extremely restricted subset of the 
physically-realizable equilibrium states of our thermally-balanced KS slab, a point not pre- 
viously known. Among the free constant parameters of this set of solutions, 70, cto and 
us define the fixed basic thermal properties of a specific fluid. Another set, B$ and Go, 
fixes the background field. The fluid is a perfect electrical conductor, so its total mass M 
frozen into the field is a constant in time. Clearly, M can be prescribed arbitrarily for a 
given set (70, a , Ks, B , G ) that specifies the nature of the fluid and the intensity of the 
background field. If the desired equilibrium state is an everywhere- analytical solution, we 
require 7 = 7 C = 2/3 which, by Equations (|20|) . determines (3 = (3 C where 



l2n>y k B 
a mBQ 



ft = ^=?. (40) 



completely fixed by the given set (70, cto, us, B , Go). By Equation f )24|) . (3 = (3 C also defines 
a unique equilibrium total mass 

M c = ^v%, (41) 

fixed independent of the above arbitrarily prescribed M. Hence, unless the prescribed M 
happens to have the value M c , the prescribed fluid cannot be in an everywhere-analytical 
equilibrium. 

To summarize, if we prescribe M = M c , M c fix by the given set (70, a , k s , B , G ), the 
fluid in equilibrium has an infinity of everywhere-analytical states available, of different Y 
dependences, depending on the central temperature #(0) = d\ as a free parameter. Each of 
these equilibrium states is the end state of an evolution from some prescribed initial state, 
governed by the time-dependent extensions of the static equations ffTBl . (fi~4D and ffT9|) . The 



-16- 



total mass M = M c is a constant in that evolution. The central temperature 6\ defines 
a minimum amount of internal energy retained for each piece of the equilibrium plasma. 
The 6>! = equilibrium state is the lowest energy state, a preferred end state in the sense 
that any higher-energy 9\ 7^ equilibrium state may evolve in response to perturbations 
by losing energy. If 9\ 7^ 0, the equilibrium D is finite, and P = D9\ — 1 at Y — 0. As 
9i — > parametrically, D(Y) as a function of Y becomes unbounded at Y = 0, diverging as 
D m Y~ v with a constant power < v < 1. In that limit, D(Y) not a Dirac delta function 
in the neighborhood of Y = 0. That is to say, the integral of D(Y) over a small interval 
— Y < Y < Y goes to zero as Y — > 0. 

The interesting new point is that, for all other prescribed M 7^ M c with /3 7^ /3 C , a 
contradiction arises if we demand that the equilibrium state is everywhere analytic. Neither 
the boundary conditions nor the freedom to arbitrarily prescribe M can be compromised. 
The physics of the governing ordinary differential equations (ODEs) are also not in question. 
Physics avoids the implied contradiction by admitting disc ontinuous solutions, kno wn as 



weak solutions to ordinary and partial differential equations (ICourant fc Hilbertlll962l ). The 
continuous parts of a weak solution are described by the governing ODEs, whereas the same 
physics governs the discontinuities by the integral version of these ODEs. The equilibrium 
states of a fluid specified by a given set (70, a , k s , B q , G ) are thus generally singular. Using 
set-theory language, the infinite set of continuous equilibrium states with M = M c is a subset 
of measure zero of all the physically admissible states. 



3. Equilibrium with a cold mass sheet 

The inviscid, compressible, supersonic flow past a fixed obstacle is a familiar example 
of inevitable singularities in solutions to partial differential equations. The flow must be 
sub-sonic in the vicinity of the obstacle to enable the fluid to self-communicate and navigate 
steadily around the obstacle. A shock front must generally be present where the supersonic 
flow changes discontinuously to a subsonic flow governed by the integral version of the ideal 
hydrodynamic equations, the Rankine-Hugoniot equations. This surface of discontinuity 
can be resolved into a continuous layer of finite thickness if we introduce additional physics 
significant only in that layer. For example, the introduction of a weak viscosity resolves 
the shock front into a layer of a thickness determined by the strength of the viscosity. The 
singular KS-slab solutions presented in this Section are analogous. The everywhere-analytical 
continuous solutions, requiring a single fixed total mass M = M c , are physically just too 
restricted. We first survey the singular solutions to understand their physical meaning. In 
the next Section, we address the additional physics one might consider to account for the 
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internal structure of the singularity. 



3.1. The complete absence of equilibrium 

Before proceeding, we should point out the case of complete absence of equilibrium 
for a KS-slab. If radiation is absent, in the limit of 7 — > 00 in Equation (jT^|) . there is no 
steady state because the unbalanced heat input must everywhere increase the temperature 
in time. In the opposite extreme, if heating is absent, with 7 — > 0, there is only one singular 
steady state. The fluid simply radiates away what heat it happens to possess to reach zero 
temperature, collapsing its entire mass M into a sheet of zero thickness; see the isothermal 
solutions in Section 2.2.1. 

When both radiation and heating are present, the total mass M critically determines 
whether energy balance is possible. In equilibrium, the hydrostatic pressure P(Y) supporting 
the fluid weight along B must be monotonically decreasing outward from the center of the 
KS slab. With the fixed heating of the form in Equation (IT51) . this stratification tends to 
produce a dominance of radiation around Y = and a dominance of heating in the far 
regions. Therefore stratification must be consistent with a temperature 9{Y) monotonically 
increasing with distance away from Y — 0, such that its inward thermal conduction brings 
the excess heat of the far regions to be radiated away by the dense slab center. 

Too little mass results in poor radiational loss, and no equilibrium is possible as the 
unbalanced heat-input keeps the temperature everywhere rising monotonically in time. Too 
much mass results in the opposite effect, the inability of heating to keep up with the ra- 
diational loss due to a high-density massive slab center. We have two possibilities in this 
case. The entire mass may collapse into single massive sheet in a time-dependent evolution, 
never attaining a steady state and leaving vacuum outside the sheet. The other possibility 
is that this collapse takes place but does not involve the entire mass. The collapse takes 
out just a fraction of the total mass so that the rest of the fluid in the outer layer achieves 
a spatially-extended steady state. These general expectations are met in the quantitative 
analysis of our KS slab given below. 



3.2. The discontinuous KS-slab solutions 

Suppose the KS slab contains a central mass sheet with a total mass m < M. The 
monotonically increasing B z = B H(Y) jumps from a value — H 2 to H 2 upon crossing Y = 
from negative to positive Y. This jump defines a discrete sheet current in the x direction 
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producing a Lorentz force that supports the weight of the sheet, 

mg = |° # 2 , (42) 

similar to formula ( I24p . The density distribution D(Y) is a Dirac delta function in the vicinity 
of Y = that accounts for m. In contrast the pressure P(Y) must be everywhere finite, 
given by Equation ffT3|) to ensure pressure balance with the magnetic field, with -P(O) = 1 
at Y — 0. The transformation — oo < K < oo — > — y/]3 < H < yT? is thus discontinuous, 
with the single point Y = mapped into the range — H 2 < H < H 2 , and < Y 2 mapped 
into (3 < H 2 . 

By converting an integration with respect to Y into its equivalent integration with 
respect to H, we need to allow for the presence of a delta function where one occurs. The 
formula ( 1371) for the total heat input % remains valid in the presence of the sheet of mass 
m. The total radiative loss needs to be re-calculated, 

K = 




:i;2 (2^ 2 + Hi). (43) 

Global energy balance 1Z = % is then expressed by 

P 2 = (3 7 -2) 1 / 3 V ^, (44) 



relating 7 to H 2 , or, m given by Equation ( 1421) . Let us fix the set (70, «o, Ks, B , Go) to 
specify a fluid of interest. Setting the eigenvalue condition 7 = 7 C yields H 2 = m = and 
M = M c , recovering the everywhere-analytic solutions. Conversely, if M 7^ M c , use Equation 
(|24|) to fix (3 in terms of the given M, which goes on to determine a 7 7^ 7 C by Equation (|20l) 
and a consistent H 2 by Equation 1)441) . The mapping —00 < Y < 00 — > — \ff3 < H < a/5 
requires H 2 to be positive and bounded above by \ffi~. Whether our M 7^ M c construction 
yields a physical solution depends on whether H 2 meets the aforementioned requirement. 

To proceed, introduce a second particular value of 7 = j' c = 1.0. This value of 7 by 
Equation ( 120]) determines a value (3 = (3' c , that defines a second critical mass M' c by Equation 
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( I24p . The nature of the equilibrium, or absence of equilibrium, depending on the total mass 
M, can now be distinguished into the ranges : M > M c , M = M c , M' c < M < M c , and, 
M < M' c . 

If M > M c and 7 < 2/3, Equation (144]) gives a negative H 2 in violation of H 2 falling in 
the range < H < y/JJ, so equilibrium does not exist even when we allow for a central mass 
sheet. There is too much mass and the heating cannot keep up with the radiational loss. So 
the entire mass must collapse to a single mass sheet located at Y = 0. 

If M = M c , we have 7 = 7 C = 2/3, H 2 = 0, and the equilibrium solution is everywhere 
analytical. 

If M' c < M < M c , 7 C < 7 < j' c , < H 2 < a/5 and the equilibrium solution exists, 
discontinuous at Y = where a m mass-sheet is located. 

If M < M' c and 7 > 7 C , then Equation ( jUJ) gives a. H 2 > \fj3, in violation of if 2 falling 
in the range < H < y/]3. No equilibrium solution exists, even if mass sheets are allowed, 
the case of the total mass M being too small for radiational loss to balance against heating. 
So the KS slab can only remain in a time-dependent state with temperature increasing with 
time everywhere. 



3.3. The M' c < M <M C Weak Solutions 

To construct an explicit singular solution in the range M' c < M < M c , apply boundary 
condition (l2~3l at H = a//3 (Y — > 00) to thermal flux T given by Equation ( 1331) : 

F(H) = ^(^-H) [H (H + v^) - P (2 - 3 7 )] , (45) 

requiring Hi = — 1/3 3//2 (2 — 37). This thermal flux vanishes at Y — > 00 but is dependent on 
7 and generally not zero at Y — 0, the point corresponding to H = H 2 . Figure 7 displays 
the functional forms of F(H) constructed, for several representative 7 values. We call these 
functional forms unmodified, having analytically extended each of them from the physical 
range of interest < H < \//3 into the negative if -range. 

The 7 = 7 C = 2/3 thermal flux is unique, with H 2 = m = and satisfying boundary 
conditions (122]) and ( J23l) at H — 0, y/~f3. So T{H) is physical and defines an infinity of 
continuous solutions in the full range — \//3 < H < y/]3, ie, — 00 < Y < 00. These solutions 
all have M = M c differing in the central temperature 6(0) = Q\. 

For 7c < 7 < 7 C , the solution is singular, obtained by first truncating off the function 
J~(H) at the limit H = H 2 of the range H 2 < H < a//3 that is now to be identified with 



-20- 



< Y < oo, for H 2 given by Equation ( )44j) . In Figure 7, the lines H = H 2 are drawn from 
the abscissa to the respective J~{H) curves in the range 7 C < 7 < j' c - The modified J~{H) 
located in the interval — \ffi < H < —H 2 , identified with —00 < Y < 0, is defined by the 
transformation J-(— H) = —F(H) as illustrated in Figure 8. These two continuous solutions 
in Y > and Y < are assured to match the cold, discrete sheet of mass m related by 
Equation (142]) to the jump of 2H 2 in H(Y) across the sheet. 

A weak solution for the 7 = 0.7 is shown in Figures 9. This solution is obtained by 
integrating Equation ( 1331) with respect to H for 9(H), using Equation lIHTl) . over the range 
H 2 < H < a//3. From it we obtain the symmetric temperature 0(Y) in < \Y\ up to a free 
central temperature that we set to zero, ie, 6(0) = Q\ = 0, where the cold sheet of mass m is 
located. However, this zero temperature is associated with a finite thermal conductive flux 
entering the sheet from both sides, in evidence in Figure 8 showing J 7 ^ at H = H 2 . In 
other words, 0(Y) goes to zero on either sides with an infinite gradient such that T is finite 
and nonzero at Y — 0. 

The central sheet should not be thought of being cold with zero temperature but as the 
limiting case of a layer cooling to zero temperature, and collapsing during an evolution to 
achieve equilibrium. The collapse proceeds with density D becoming a Dirac delta function 
to account for the mass m and with P holding steady at unity. The right hand side of 
energy-balance equation ffT9|) integrated over the thickness of the sheet gives a constant in 
time, the steady net loss of energy by the collapsing sheet balanced by the steady inward 
thermal fluxes on both sides of the sheet. During such a collapse, D increases in the sheet 
to become a Dirac delta function, temperature decreases to keep P = 1, producing a steady 
net loss. This is the nature of the 2/3 < 7 < 1 equilibrium states, in each case taking the 
macroscopic physics of Equation ( Tl9|) to remain valid for all time in the interior of the mass 
sheet. In the next Section an example is given for a collapsing sheet that has a distinctly 
different energy process from the exterior continuous parts of the SK-slab. 



4. The interior physics of a discrete mass sheet 

To complete our study, we present two alternative interior-models for a mass sheet 
matching an exterior continuous solution of our KS slab. As a sheet collapses with decreasing 
temperature and unbounded growth in density at a fixed total mass m, we expect the 
optically-thin assumption to break down. The nature of heating would also be changed. 
Heating in the prominence is poorly known or understood, which is the reason for adopting 
the ad hoc form given by Equation ffT8|) . used in earlier prominence modeling works. 
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For illustration purpose we consider the case of a finite-thickness mass sheet whose 
interior is optically thick to all radiation, subject to no heating, but maintained by thermal 
conduction to be isothermal, allowing for the uniform temperature to be as low as it happens 
to be at the interface with the exterior region. In the exterior, described by our KS-slab 
solutions, the total heat generated is exactly balanced by the total radiative loss. No thermal 
flux enters mass sheet. When this sheet is at a sufficiently low temperature, the degree of 
ionization would be so low that the frozen-in condition must break down. Then, resistive 
flows of weakly ionized fluid across the "supporting" magnetic field can take place, the second 
alternative interior-model that takes us a step closer to the observed prominence with its 
descending threads and rising bubbles. 



4.1. Composite KS-slabs 

Consider the case 7 = 9/16 in Figure 7, as an example of M c < M, i.e., 7 < 2/3. The 
unmodified thermal flux J-{H) has two zeros in the range < H < y/]3, one at H = yf]3, 
imposed by the boundary condition at Y — > 00 and the other at a point at < H — if 3 < \f^. 
This means that the heat generated in the outer layer H3 < H < \fj3 is entirely radiated 
away in that layer. In other words, there is enough mass in the system to create an outer 
layer that is by itself in global energy balance. This layer provides no energy to be transferred 
by thermal conduction to the inner layer, whatever the latter may be. Integrating Equation 
(|45|) we obtain 0(H) up to a free constant 9(H 3 ) = #1, the temperature at the interfaces 
H = ±H 3 . In a time-dependent cooling of a M c < M slab, depending on the nature of 
its initial state, it is possible for such an outer layer to first form and and then stabilize, 
symmetrically placed on the two sides of Y = 0, filling up the regions \H\ > if 3. The rest of 
the mass m in \H\ < if 3 may then exist in two forms which we now specify. 

If we define the mass m as a collapsed, discrete mass-sheet at zero temperature, we have 
nig = -^Hz in analogy to Equation fj4*2l) . Then, the monotonically increasing B z = B H(Y) 
jumps from — H 3 to H 3 across Y — 0, with the point Y = mapped into the interval 
—if 3 < H < H3 and the intervals < Y 2 into the iintervals iff < H 2 . We must set the 
boundary condition 6(±Hs) = Q\ = for the external solution. A static solution of this kind 
with 7 = 9/16 is displayed in Figure 10. The energy equation (TT9l) applies to the exterior 
solution but not to the interior of the sheet. 

The mass m can also exist as a symmetric, isothermal slab of a prescribed temperature 
6*i, of a finite width occupying — Y3 < Y < Y3, matched continuously into a pair of external 
solutions of our KS-slab occupying symmetrically in —00 < Y < — Y3 and Y3 < Y < 00 
subject to the boundary conditions that the exterior temperature 6 = 61 at Y = ±Yjj. First 
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prescribe 6\ and pick up the value of H% from Figure 8 for 7 = 9/16. The exterior temperature 
9(H) in H 2 < H 2 < (3 is then determined. The picked value of H 3 determines m, sufficient 
for constructing an isothermal solution that locates H = ±H 3 at a pair of points Y = ±Y 3 
in physical space. This isothermal solution then self-consistently matches the outer solutions 
at Y = ±Y"3 with 0(±H^) = 0\ and J-(±H$) = 0. The construction is then complete, giving 
continuous distributions P(Y), D(Y), 0(Y) and H(Y) over —00 < Y < 00. Figures 11, 12, 
13, 14 display 4 solutions computed in this way, all sharing the same M and m but for a set 
of progressively smaller temperatures 9\ as indicated. Although a dynamical time-dependent 
theory is outside the scope of our study, our examples mimic the cooling of the central core 
of a fixed mass m, holding the total mass M also fixed. As the central isothermal core cools 
quasi-steadily, that core collapses into a slab of diminishing width, approaching a Dirac delta- 
function in the limit of 61 — 1 0. Such a collapse proceeds first out of its optically-thin state 
into an optically-thick state, and then to extremely low temperatures when the assumptions 
of high ionization and the frozen-in condition must fail. 



4.2. Resistive vertical cross-field flows 

In an extremely cold and dense fluid, the electrical conductivity cannot be taken as 
infinite by two effects, the low degree of ionization and the narrowing of collapsing layer 
to such a small size that the weak but finite electrical resistivity has become important. If 
resistivity is significant, the Kippenhahan-Schluter slab can be extended to include a steady 
vertical resistive flow. Consider the single-fluid resistive hydromagnetic equations ([I]) and 
02]), taking the resistivity 77 uniform in space for simplicity. In the presence of the simple 
steady vertical velocity 

v=[0,0,v(j/)], (46) 

the rate of change of momentum on the left hand side of Equation ([I]) vanishes identically 
and we recover the force-balance equation ([3]). We also recover the steady energy transport 
equation (J5J). The mass conservation is automatically satisfied. Therefore, each of the mag- 
netostatic solutions investigated in this paper is consistent with the above steady, vertical, 
resistive flow with 77 7^ 0. 

The induction equation (J2J) reduces to a simple second-order ODE which can be inte- 
grated once to give 

^ ] --iw + ^ (47) 

where vq is an integration constant and Lq is the density scaleheight. In all our solutions 
tS — > as Y — > 00. Therefore, we set vq = for a system static in the far regions. Combined 



-23- 



with magnetostatic equation (JHJ), we have the velocity 

v(y) = -£d. (48) 

In this resistive flow the energy released from current dissipation is conserved by the release 
of gravitational potential energy associated with the steady fall of matter. The prominence- 
drainage event of Liu, Berger & Low (2012) is thus interesting. The drainage liberates 
gravitational potential energy. If that drainage is across field lines, the liberated energy 
manifests, not as an increase in kinetic energy, but as dissipative energy that fuels the 
radiative loss of the prominence plasma. 

Suppose rj is small but not zero. Then to a first approximation we neglect rj and v. The 
results of Section 3 have the following significance. Suppose M > M c , static equilibrium is 
achieved in a relaxation by withdrawing a specific portion m of the given total mass M into 
a collapsed central sheet. This process implies that the density D grows unboundedly at 
Y = 0. At some point, for any given rj small but not zero, the velocity v becomes significant 
at Y = leading to a breakdown of the frozen-in condition in spite of the assumed smallness 
of r\. 

For example, consider this resistive flow taking the form 

"M = -m^w/' (49) 

with D given by Equation f l26|) at the uniform temperature T = 9\Tq. We have a vertical 
down flow that is maximum at Y = but decreasing with \Y\ exponentially. Use the 
normalization T = 10 4 K so that L « 300 km. The Spitzer electrical conductivity a = 
2 x 10 7 T 3 ' 2 s^ 1 gives a resistivity rj m 3 x lO 6 ^ . The peak resistive flow at Y = is 
then of the order of nl .„ ^5x 10 _2 ^o. To produce a flow of the order of 1 kms~ l , would 

require (3 large and Q\ small, e.g., /3 ~ 10 3 , Q\ ~ 10 -2 , implying a highly sheared field and 
low temperature of the order of T oc 10 2 i^. Such an extreme condition is the reason for a 
long-standing pessimistic perception that steady resistive slippage of ionized plasmas as the 
result of the Spitzer conductivity is unimportant. The spontaneous formation of a mass- 
sheet singularity is a novel example of how precisely this kind of extreme circumstances may 
develop in the corona under which resistive slippage does become important. 

This effect is more interesting as a recurrent time- dependent process. If we neglect the 
weak resistivity altogether, the frozen-in condition demands for an equilibrium state with 
a singularity. This singularity generates a resistive flow of an unbounded magnitude until 
the singularity formation is stopped by current dissipation as a heat source for the cold 
fluid. The heating can increase the degree of ionization, restoring the frozen-in condition, 
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whereupon a collapse into a mass sheet is once again initiated. This time- dependent process 
is worthy of further investigation. 



4.3. Ambipolar diffusion in the KS slab 



To address the degree of ionization, a two-fluid description is needed ([Gilbert. Hansteen fc Holzer 



20021 ; iGilbert. Kilper fc Alexanderll2007l ). Such a description is necessarily 2D, quite outside 



the scope of our study, but the following is an interesting closing point to make. 

In a ID slab of weakly ionized plasma, the dominant neutrals are kept from slipping 
freely across the field by collisions with the ions that tend to be tied to the field. The 
system cannot remain ID because there will be diffusion of the neutrals in both Y and Z 
directions. If we concentrate on the vertical force balance in the central region of a KS- 
like slab, the following order of magnitude consideration adapted from Spitzer (1978, page 
294, eqn. (13-55)) is instructive. Take the ions to be rigidly tied to the locally horizontal 
field as the neutrals fall steadily vertically across the field lines in the presence of a mutual 
frictional force. For the neutrals as a fluid, its weight is balanced by the upward frictional 
force acting on it, denoted by F z ^ n . The ions as a fluid experiences the opposite frictional 
force F z ^i = —F z ^ n as well as the upward Lorentz force that effectively supports the total 
weight. Use the Spitzer approximation, 

F Z)i) n = -F z ^i = pip n — v(y), (50) 

where v > is the speed of the neutral fluid falling relative to a stationary ion fluid with 
densities p n and p i7 respectively. The parameter c is a thermodynamic constant not greatly 
sensitive to temperature, ranging in value from 1CT 10 to 10~ 9 for a hydrogen gas at temper- 
atures in the range 10 — 10 3 f^, and, we recall, mo is our notation for the mean molecular 
weight of the fluid which we shall take to be the hydrogen mass. We also have the vertical 
balance of forces for the two fluid components: 

p n g = F Zjij n, (51) 

Pig = F z>n>i + Bl — . 52 

ay 

Summing these two equations gives the one-fluid Equation (IT41) following the canceling of 
equal and opposite frictional forces. 

Substituting for the frictional force in Equation ( 15"T|) given by Equation ( 1501 . we obtain, 
after a canceling of p n as a common multiplicative factor on the two sides of the equation, 



the speed of fall is independent of p n : 
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•<»> = "S- (53) 



entirely determined by the number density N{ of the ionized fluid. The speed of fall depends 
on gravity but it is steady and not of the magnitude of a free-fall speed. A slippage velocity 
of 2 km s~ l is obtained if the ion number density drops to below 10 5 cm~ 3 . Again, this speed 
of fall would not be significant if not for the principal result of the paper. If no such poorly 
ionized fluid form, the frozen-in condition would obtain. Then, if the mass frozen in the the 
field is greater than M' c , the collapse of the core of the KS slab would bring the temperature 
at the center to as low as it takes for the number density of the ions to be reduced to the 
level where the frozen-in condition is invalid. Heating is produced by the friction force and 
the flow can be thus quenched when the ionization level is increased. What is interesting 
is that if the frozen-in condition returns, and the total mass still exceeds M' c , cooling and 
collapse of the optically thick core recur. The condensation is thus recurrent and not static 
for these theoretical reasons. 



5. Discussion 

The vertical descent of prominence threads at less than free-fall speeds has been inter- 
preted in several ways. The large-scale twisted field in the prominence cavity may be evolving 
with flux-tubes falling bodily together with their frozen-in plasmas as t he result of intermit 



tent magnetic reconnections occ urring somewhere else along a flux-tube ( ILow &: Petridl2005 



Petrie fc Lowll2005l ; I Chad |2 01 Of ), van Ballegooijen and Cranmer (2010) proposed that the 
prominence plasma is flowing under the frozen-in condition in a 3D tangled global field that 
has a strong vertically-oriented component. In contrast, here we propose a spontaneous 
breakdown of the frozen-in condition to allow a resistive flow across a local horizontally- 
oriented field. What seems clear is that all these dynamical processes may occur simultane- 
ously. 

Our KS-slab model shows that the nonlinear coupling of force-balance with energy- 
balance under a rigorous imposition of the frozen-in condition must generally produce a 
cold, zero-thickness sheet of mass if static equilibrium is to be attained. This mass sheet 
contains a discrete electric current that generates the Lorentz force supporting the weight of 
the sheet. In a real plasma of a high but finite electrical conductivity, this discrete current 
must dissipate resistively. In other words, the typically high degree of the frozen-in condition 
must physically break down within the mass sheet. In a time-dependent relaxation to such 
an equilibrium, a mass sheet would first develop during the time when the frozen-in condition 
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applies approximately. Then, in that developing mass sheet, the frozen- in condition soon 
breaks down to produce a downward resistive flow across the supporting field. Yet, the 
frozen-in condition is not broken down permanently. The implied resistive heating in the 
sheet, fueled by gravitational potential energy released by the falling, poorly- ionized plasma, 
can bring back the frozen-in condition, but only to recreate the circumstance under which a 
thin, cold sheet is again inevitable. This recurrent spontaneous breakdowns of the frozen-in 
condition is a physically attractive hypothesis to explain the constantly restless interior of a 
quiescent prominence as revealed by SDO/AIA and Hinode observations. 

Gravity, anisotropic Lorentz force, optically-thin radiation, heating, and anisotropic 
thermal conduction all play essential roles in creating what seems quite general circumstances 
under which thin collapsed mass sheets must form. The particular optically-thin radiative 
loss and heating described by Equations ( 1T7|) and ([IS]) are artificial, but they capture basic 
properties sufficiently well for us to use the KS-slab to demonstrate the inevitability of mass- 
sheet singularity with mathematical rigor. Our demonstration is surprisingly complete. The 
mass-sheet singularity in our KS-slab is resolvable into a sheet of finite width using an exact 
solution for a cross-field resistive flow to account for the breakdown of the frozen-in condi- 
tion. Both the static and time-dependent problems of this complex coupling are formidable, 
especially for 2D and 3D systems. Direct numerical simulation holds the only hope for mak- 
ing progress but to be fruitful this approach must start with a physics-based hypothesis like 
the one presented here. Much further work is needed to treat physically more realistic en 



ergy processes for applications dire ctly to the prominence (JGilbert. Hansteen fc Holzerll2002 



medium ( JMac Low et al.lll995t 



Gilbert. Kilper fc Alexander! 120071) and other astroph ysical systems, notably the interstellar 



Parrish & Stond 120051 ). 



Observation remains the s p earhead in the inv e stigation of quiescen t prominences (jBerger et al 



2008 



2010. 



2011 



0; 



0; 



Chad 120101; iGaizauskasI Il998l: iHeinzel et al.l 120081; lLabrosse et al. 



2010; 



! 994; Martin fc Panasencol 



Lites et al.l 



201 



Mackav et al. 



2010 |; lTandberg-Hanssenlll995l: 



20K 



Labrosse et al 



201 



0; 



201 



Martim_Bi 



Lites et al. 



Mackay et al. 



2010: 



2010 



imoria fc P. W. Tracada; 



Okamoto et al 



200 



van Ballegooijen 



7. 



2008 



Okamoto. Tsuneta fc Berger 



20041 ). Presently not all 



the observed properties in prominences can be theoretically understood or understood in 
a unified way For example, the sub-arcsecond threads, revealed by ground-based, high- 
resolution, static H a images of on-disk prominences, have remained difficult for theory to 
unify with the space observations of prominences at the solar l imb from Hinode and SDO/AIA 
rtzirker. Engvold & Marth]|l998l ; Em". Engvold & Wiikil2003[ Ein et al-lbooi ISu et alJboilh . 
These two views of prominence structures have also not bee n reconciled with the polarimet 
ric observations of prominence and corona 
Casini fc Judgdll999J [Judge. Low fc C 



2008 



Truiillo Bueno et al.ll200a iDove et al.l 1201 11 ) 



vasim 



magnetic fields (ILerov 



2006 



1989 



jQpez Ariste fc Casinill2003l ; iTomczyk et al 



Casini et al. 



20031: 
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Observations can pose good physical questions to be addressed theoretically in their 
own rights, without necessarily aiming at a comprehensive explanation of all the rich obser- 
vations. In the face of many unanswered questions, our series of paper seeks to build up our 
intuition for seeing the complex observations in terms of the elementary physical principles. 
Previous works of this kind have addressed magnetic support of the prominence weight under 
the frozen-in condition using 2D and 3D fi eld models with energy transport greatly simpli- 



fied or ignored for the sake of tractability (Aulanier. DeVore fc Antiochos 



Manchester et al 



2004 



2004 



van Ballegooijenl I2004J ; IWu k Low 



2002 



Chae 



Low k HundhauserJfl995l;lLow k ZhanelboollLow k Petridhood: Lopez Ariste et al. 



Petrie k Lowlbood; IPriest. Hood k Anzerl Il989l ; iRegnier k Amari 



2010; 



2006; 



19871 ). In complement, models of ID sophisti- 



cated energy transport along field lines, prescribed to be fixed in space, made progress 
with insights into the condensa tion processes by suppressing the dy n amical interaction of 
the field with energy transport (IKarpen. Antiochosfc Klimchukl 120061 ; iKarpen k Antiochos 



2008 



Luna. Karpen k Devorell2012l ). Heasley and Mihalas (1976) pioneered studying the 



coupling of radiative transfer to the equilibrium forces in a slab model. Allowing for gravity 
and a heuristic approximation of the Lorentz force, this interesting work encountered an 
inevitable presence of collapsed mass-sheets, perhaps a rudimentary form of the singularity 
we found in our KS-slab model; see also Anzer k Heinzel (1999). 

Collapsed mass sheet may be a common occurrence in the corona. All it takes is a failure 
to bring energy to where it is needed whenever a particular mode of energy transport becomes 
over constrained by the coupled requirement of balancing the forces in the gravitationally 
stratified fluid. Let us conclude with the following intuitive physical picture that has emerged 
from our study. 

In our KS-slab, there are two thermal-insulation mechanisms that can keep the promi- 
nence cool against the million-degree hot, thermally highly-conducting corona. The magnetic 
field is an insulator against thermal conduction across it. Not so well recognized is the self- 
insulation of the prominence plasma along a flux tube. In a sheared U-shaped flux tube, 
the frozen-in mass is gravitationally stratified along the tube so as to progressively radiate 
away the two thermal fluxes conducted inward from the high-temperature upper parts of the 
tube in the corona. If the total mass is large enough, the steady state may be one in which 
these thermal fluxes are depleted before the bottom of the tube is reached. Then the bottom 
contracts into a cold dense mass simply from want of heat. The degree of ionization of such 
a core would fall below a threshold whereupon the frozen-in condition must break down. 

The prominence/corona transition is similar to the corona/chromosphere transition of 
the solar atmosphere but that analogy is just imperfect. In the solar atmosphere, the bottom 
of the chromosphere never cools below about 5000°/^ because the photosphere beneath it 
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has no bottom and is heated actively by the solar luminosity coming through from below. 
If the prominence plasma sitting at the bottom of an U-shaped tube has no internal heat 
source and is self-insulating, it would simply cool to an equivalent of the solar photosphere 
with a very different property, one with such high density and low temperature as to break 
the frozen-in condition. This massive core would then flow resistively across the field into the 
adjacent flux tubes. No steady flow is expected, just a recurrent process of the breakdown 
alternating with restoration of the frozen-in condition by resistive heating. 

In realistic 3D situations, the super-cooled core of a prominence plasma slips resistively 
out of its optically-thin chromospheric cladding into an adjacent flux tube. There it can 
either continue to slip resistively into other tubes further away or else become so ionized in 
the new hot environment as to be frozen into the new tube it has just occupied. A sufficiently 
massive super-cooled core can thus fall through the bottoms of successive U-tubes leaving 
behind in each of these tubes ionized plasmas at chromospheric temperatures, until the 
former runs out of mass to remain a collapsed core. A vertical thread of chromospheric 
plasma is naturally produced. The movie of such a process made at the cool chromospheric 
temperature would then show such a thread forming from out of nowhere with a general 
downward descent. The presence of U-shaped flux tubes is the pre-requisite of this process. 

Although such massive, poorly ionized cores are time-dependent objects, this picture ar- 
gues for the presence of prominence plasmas that are more massive, more dense, much colder, 
and less ionized than currently being contemplated. Perhaps the canonical observation- 
based numbers describing the prominence, electron density of 10 10-12 cm~ 3 , temperature of 
~ 5000°K, degree of ionization 0.1 — 10%, are only the properties of the optically-thin, chro- 
mospheric cladding of the not yet observed, transient, optically-thick cores. The tangled-field 
models of van Ballegooijen and Cranmer (2010) readily give peak densities well in excess of 
the aforementioned canonical values. These models assume a moderate field of only 10G, 
whereas even larger densities are implied for fields as high as 60G given by polarimetric 



observations ( ICasini et al.l l2003t lLopez Ariste &: Casinil 120031 ). Our KS-slab suggests that 



even more extreme densities might occur commonly in prominences. 

We emphasize that upward and downward motions accompany this process. The un- 
loading of a flux tube to give mass to the one below results in the former rising by its resultant 
overcompensating Lorentz force. Correspondingly, the mass-receiving tube below sinks heav- 
ily until the Lorentz force of its stretched field can hold it swaying into position. Here is an 
instance of a flux-tube of fluid rising or falling bodily, carrying its frozen-in magnetic flux 
along with it. 

Our hypothesis suggests that the cavity magnetic field supporting the prominence is 
"porous" in spite of the high electrical conductivity. The observations of Liu, Berger & Low 
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(2011) thus have the significance as a quantitative statement of how porous that magnetic 
field can be, namely, that an order of magnitude more mass than retained by the prominence 
at any one time can pass resistively and agitatedly down the cavity field in a day. 

In t he second paper we will get a clearer picture of the role of the magnetic field in this 
process (ILow et al.ll2012l ). The field actually is the dominant dynamical player, for it can, on 



its own, create and restore the frozen-in condi tion in a permanent restlessness (jParkerlll994 



Janse fc Lowil2010i ; |Janse. Low fc Parkerll2010l ). This mechanics may be basic to the proposal 
that the long-lived macroscopic prominence is a form of magneto-thermal convection in the 
corona proposed by Berger et al. (2011). 

We thank Yuhong Fan, Michael Goodman, Phil Judge and Bruce Lites for discussions. 
The National Center Atmospheric Research is sponsored by the National Science Foundation. 
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Fig. 1. — A Hinode/SOT image in Call H emission of the fine-scale vertical threads in a prominence 
observed on October 3, 2007 at the solar limb. 




Fig. 2. — White-light total eclipse of March 18, 1988 showing a well-formed, bright coronal helmet- 
streamer extending from the north-east limb, with a localized, low-density cavity at the helmet base. 
Within that cavity a spatially-unresolved, quiescent prominence appears as a bright blot. 
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Fig. 3. — The isothermal Kippenhahn-Schliiter prominence slab displayed with graphs of density 
D(Y); vertical field component H(Y); and, the shape of the bowed field line projected on the Y — Z 
plane described by Z(Y), as described in the text. 
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Fig. 4. — Thermally-balanced continuous solution with /3 = 1.0, Go = 0, showing the sharply 
spiked density D(Y); rapid rise of temperature 6{Y) from a low value at the origin; the vertical 
field component H(Y); and the shallowly bowed field line represented by Z(Y) in the Y — Z plane. 
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Fig. 5. — Thermally-balanced continuous solution with (5 = 5.0, Go 
Figure 4. 



0, in the same format as in 
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Fig. 6. — Thermally-balanced continuous solution with /3 
as in Figure 4. 



130.0, Gq = 5.0, in the same format 
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Fig. 7. — Unmodified thermal flux J-(H) for different values of 7 indicated over the range —y/]3 < 
H < yffi. Except for the antisymmetric 7 = 2/3 graph, vertical lines H = Hi are drawn from the 
abscissa to the respective curves as described in the text. 
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Fig. 8. — Modified thermal flux J-{H) for different values of 7 indicated over the range —\f]5 < 
H < \f~fi. Except for the antisymmetric 7 = 2/3 graph, the modified F{H) in H > is defined 
only for H > H2 truncating off the —y/]3 < H < H2 portion of the unmodified J-(H) shown in 
Figure 7. The modified J-(H) in H < is then defined for —y/]3 < H < —H2 by an antisymmetric 
reflection of the part so constructed for H > H2 ■ 
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Fig. 9. — The density D(Y) with a central Dirac delta function and temperature 6(Y) of a 7 = 0.7 
thermally-balanced solution. The delta function describes a collapsed mass sheet corresponding to 
the jump of 2H2 across the origin of the vertical field component HiY). The thermal flux T for 
7 = 0.7 shown in Figure 8 indicates that 9(Y) — > with an infinite gradient as Y — > 0, delivering 
a non-zero thermal flux into the mass sheet as described in the text. The constant #0 is only a 
normalization to fit D(Y) and 9{Y) into the same graph. 
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Fig. 10. — A 7 = 9/16 thermally-balanced solution with a collapsed mass sheet in the same format 
as in Figure 9 except that, not resolved at the scales of the figure, 6 — > with zero thermal flux 
at Y = as discussed in the text. In the continuous part of the solution, the total heat input is 
balanced exactly by the total radiative loss. 



-42- 



1,0 



Two— component solution; gamma=9/16, beta=1.0, G_0=10.0 
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Fig. 11. — The distributions of the vertical field component H(Y) of the two fluid-component 
thermally-balanced solutions displayed in Figures 12, 13, 14, distinguished by the central temper- 
atures 9\ of these solutions as indicated, measured in a common unit. Each solution is everywhere 
continuous, having the same total mass M and thus the same range —yf$ < H(Y) < \J]3, with 
an isothermal solution of temperature 9\ in the range — H$ < H < #3 and a 7 = 9/16 thermally- 
balanced exterior solution in H3 < \H\. The parameter H3 is the same constant for all three 
solutions, as described in the text. 
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Two— component solution; 7=9/16, ,3=1.0, G„=10,0 



0.0 



, , I I I I I I I 


, I I I I I I I I 
\ d e, 

y s/io \ 









-10 




Y 



10 



Fig. 12. — The broad distributions of density D(Y) and temperature 0(Y) with central temperature 
9\ = 4.0. The latter shows clearly the location, —H3 < H < H3, of the inner fluid component by its 
spatially uniform temperature that continues into a 7 = 9/16 thermally-balanced exterior solution 
located in H3 < \H\. 
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Fig. 13. — The distributions of density D(Y) and temperature 0(Y) with central temperature 
9\ = 2.0 showing a narrowed isothermal inner fluid component, in the same format as in Figure 12. 



-45- 



Two— component solution; 7=9/16, ,3=1.0, G„=10,0 




Fig. 14. — The distributions of density D(Y) and temperature 9(Y) with central temperature 
9\ = 0.3 showing a greatly narrowed isothermal inner fluid component, in the same format as in 
Figure 12. The temperature continues from its isothermal part into its 7 = 9/16 with zero gradient, 
a feature not resolved at the scales of the figure. 



